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Interrogating individual molecules within bio-membranes is key to deepening our understanding of 
biological processes essential for life. Using Raman spectroscopy to map molecular vibrations is ideal to 
non-destructively 'fingerprint' biomolecules for dynamic information on their molecular structure, 
composition and conformation. Such tag- free tracking of molecules within lipid bio-membranes can 
directly connect structure and function. In this paper, stable co-assembly with gold nano-components in a 
'nanoparticle-on-mirror' geometry strongly enhances the local optical field and reduces the volume probed 
to a few nm^, enabling repeated measurements for many tens of minutes on the same molecules. The intense 
gap plasmons are assembled around model bio-membranes providing molecular identification of the 
diffusing lipids. Our experiments clearly evidence measurement of individual lipids flexing through telltale 
rapid correlated vibrational shifts and intensity fluctuations in the Raman spectrum. These track molecules 
that undergo bending and conformational changes within the probe volume, through their interactions with 
the environment. This technique allows for in situ high-speed single-molecule investigations of the 
molecules embedded within lipid bio-membranes. It thus offers a new way to investigate the hidden 
dynamics of cell membranes important to a myriad of life processes. 

Lipid membranes that constitute the surface of cells perform vital biological functions, from molecular 
trafficking to signaling^'^. Composed predominantly of organic lipids, these bio-membranes have become 
an intense area of research across many disciplines as efforts are made to understand the complex physio - 
chemical dynamics of the lipids^'^. Moreover, the cell membrane is built up of a complex mixture of phospho- and 
glyco -lipids as well as numerous proteins that interact with each other in sophisticated cascades that sustain 
biological pathways. To probe and identify non-destructively at the single molecule level the interactions on these 
membranes would be an invaluable investigative and diagnostic tool. Not only would it enable greater under- 
standing of the interaction of drugs and viruses that incorporate within the cell membrane, but also important 
fundamental life processes such as membrane budding and fusion. Furthermore many degenerative diseases, 
such as Alzheimer's, Parkinson's, cystic fibrosis and muscular dystrophy are believed to result from damage to the 
cell membrane. Probing single molecule behaviour in lipid bilayers can thus track in depth the development of 
these diseases. 

To observe in situ the dynamic changes in conformation and orientation of single lipids composing the bio- 
membrane is extremely challenging. Raman vibrational spectroscopy is ideally suited to provide information 
sensitive to the molecular species, orientation of constituent dipole moments, chemical state and local stresses in 
the form of a unique vibrational 'fingerprint'^'^. The challenge for single-molecule spectroscopy of biological 
systems is to provide a localized enhancement of the Raman signal sufficient to allow measurement with high- 
enough signal to noise, and to do so on the sub -second time scales necessary for following membrane dynamics. It 
is also important that this probe does not damage the delicate biological sample, allowing continued observations. 

Placing gold nanoparticles extremely close to a gold mirror (Fig. la) forms a construct well-suited to tightly 
confine light within a gap of nanoscale dimensions, in a 'gap plasmon'. Such gap plasmons have been used for 
sensitive spectroscopy of graphene monolayers^ and gap-localized chemical reactions^ with concentrations in the 
zeptomolar range^"^\ This gold nanoparticle-on-mirror (NPoM) system is ideally suited for enhanced vibrational 
spectroscopy due to the biological compatibility (inertness) of gold and the ability to create field volumes only a 
few molecules across^ '*'^^. Its advantage over previous work utilising metal tip-enhanced Raman (TERS) is not 
only its greater stability but also its ease of use while having a similar or even higher spatial resolution and field 
enhancement^^"^'*. 
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Figure 1 | Formation and characterization of gap plasmon sensor, (a), Schematic of gold nanoparticle (diameter <i = 80 nm) deposited upon hpid- 
alkanethiol hybrid bilayer on a planar gold surface, in solution. Lipids are POPCiDOTAP and SAM is octadecanethiol. (b), Wide-field scattering image of 
gold nanoparticles on lipid hybrid layer, (c) , Reproducible scattering spectra of the gap plasmon resonance from individual gold particles (shown in inset) 
on the hybrid lipid layer. 



In this work, as a model for the biomembrane we use a mixture of 
the lipids l-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) 
and l,2-dioleoyl-3-trimethylammonium-propane (DOTAP), in an 
80:20 concentration ratio. Phospholipids with the choUne group 
are chosen as our model because they are one of the most abundant 
lipid species in eukaryotic cells^^. This lipid leaflet is supported on a 
self- assembled monolayer (SAM) of long tail alkanethiols (octadeca- 
nethiol, ODT) to form a hybrid bilayer leaflet (HBL, Fig. la), a model 
analogue to the cell membrane^^. The HBL is kept hydrated in buffer 
solution to preserve the uniform membrane structure. At room and 
biological temperatures the POPC:DOTAP leaflet is in the fluid 
phase. The low steric forces of this phase allow for a lower mean tail 
density, greater tail disorder and high mobility of the lipid^^. This 
easier lipid movement is then investigated via the gap-plasmon- 
enhanced Raman signature. 

The gap plasmon cavity is formed between a flat Au surface and 
80 nm diameter gold nanoparticles deposited in situ on top of the 
supported HBL. Illuminating gold particles with a high numerical 
aperture objective excites the dipole mode orthogonal to the surface, 
and this dipole strongly interacts with its image in the planar gold 
surface. The resonance wavelength of this hybridized gap mode is 
exquisitely sensitive to the separation of the particle from the plane 
surface. Viewed under dark-field scattering microscopy (Fig. lb,c, 
Fig.Sl), the gap mode wavelength of the NPoM HBL is observed at a 
well-defined wavelength Ic^ 660 nm, permitting efficient Raman 
laser excitation with X^=633 nm. This gap mode wavelength 
matches our calculations for plasmonic coupling across a gap of 
5 nm, corresponding to the height of the HBL of 4.4 ± 0.6 nm 
measured by atomic force microscopy (AFM, Fig.S3), with a dielec- 
tric constant £^ ~ 2^^. The reproducibility of the gap mode is attrib- 
uted to the ordering of the HBL spacer. The variations in particle 
colour seen in Figure 2(b) are due to differences in height across the 
HBL, revealed by AFM (see Supp. Info.). The smoothness of the gold 
surface is measured to be <0.5 nm. Each nanoparticle initially dif- 
fuses around before settling on the lipid surface and becoming red 
under scattering microscopy (Fig. Ic, 1-3, Fig.Sl). Some nanoparti- 
cles however remain further spaced (green. Fig Ic, 4). We perform all 
spectroscopy on particles observed to be stationary and with the 
660 nm gap plasmon resonance. 

The stability and reproducibility of the gap plasmon field demon- 
strated above enables reproducible measurement of the POPC:DOTAP 



membrane, shown in averages of repeated time series (Fig. 2a). The 
field enhancement of the gap generates sharp vibrational modes with 
high signal-to-noise, which can be assigned to the known tail alkane 
vibrations and also the choline head group for the membrane lipids. 
A distinct feature of HBL spectroscopy from the NPoM gap plasmon 
is the high density of vibrational modes observed, attributed to the 
tight confinement of the cavity field around the probe molecules. 
Similar traits are also observed in single molecule SERS using gap 
plasmons between two nanoparticles where the field is similarly con- 
fined^^"^\ Here the field is further confined within the HBL due its 
dielectric contrast with the surrounding water. For plasmon fields 
with lower spatial confinement, such as with TERS or plasmonic 
substrates such as Klarite™, this spectral density is not observed 
(Fig.S4 and [22-25]). The rich vibrational spectrum from the gap 
plasmon construct is due to escaping the typical averaging over 
molecular dynamics in SERS, as discussed below. It is also partly 
attributed to modified Raman selection rules induced by the high 
degree of alignment of lipids normal to the planar Au surface along 
the very strong gradient in optical field^^'^^. This abundance of vibra- 
tional modes, each sensitive to the highly-polarized gap field orien- 
tated perpendicular to the planar surface, clearly enables the SERS 
spectroscopy to report minute changes in the orientation, conforma- 
tion and composition of the probed molecules^^. 

The discrimination in lipid species afforded by the NPoM gap 
plasmon is demonstrated by comparing hydrogenated DOTAP 
(20%) and deuterated POPC (80%) in the HBL (Fig. 2b). These 
time-ensemble vibrational spectra reveal a series of bands in the D- 
POPC HBL spectrum (highlighted) downshifted in frequency from 
those of H-POPC. The Raman bands at 3000 cm"^ and 900 cm"^ 
originate respectively from CH2/CH3 stretches and CH2 twisting/ 
wagging of the alkyl tails^^'^^. In the deuterated analogue, the hydro- 
carbon CD groups have a greater mass and thus a lower vibrational 
frequency that is clearly observed in the enhanced Raman spectrum. 
In addition, the D-POPC lipid species possess some vibrational 
modes unaffected by the deuterium substitution. Accordingly these 
vibrations occur at the same frequency in both spectra (C-C stretch 
1010 cm"^ and N-H stretch 720 cm"^), confirming the discriminat- 
ory ability of the NPoM SERS'°. 

The high field enhancement afforded by the gap mode allows for 
fast acquisition (<ls) of vibrational spectra so that time -dependent 
dynamics of the fluid bilayer may now be tracked. As an example of 
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Figure 2 | Lipidomics of lipid layers with gap plasmon SERS. (a), SERS spectra (1 s acquisition time, averaged over 300 subsequent spectra) for three 
separate particles on a POPCiDOTAP 80:20 hybrid bilayer. (b), SERS spectrum of D-POPC:H-DOTAP (10 s acquisition time) highlighting Raman 
modes shifted due to deuteration. (c), Selected spectra from repeated measurement of a three component lipid layer including 5% PTPC possessing 
unique acetylene signature mode highlighted in molecule. Inset: SERS intensity of PTPC acetylene mode with time (normalized to CH2 stretch intensity 
and acquisition time 1 s). All spectra in a-c acquired with 1 mW excitation at X^=633 nm. 



this (Fig. 2c) we measure the gap plasmon SERS from an HBL con- 
taining also a third lipid in low concentration (PTPC, 5%). This 
PTPC lipid contains a triple carbon bond in one tail with a distinct 
unambiguous Raman vibration at 2200 cm~\ By continuously mea- 
suring the SERS spectrum from a single NPoM over time, the dif- 
fusion dynamics of PTPC within the membrane is observed. The 
acetylene signature vibration of PTPC fleetingly appears as the lipid 
molecules diffuse into and out of the gap plasmon mode volume 
intersecting the membrane. The evolution of this signature mode 
(Fig. 2c inset) over 300 s (normalized to the 1380 cm"^ CH2 ref- 
erence) implies that single molecules can be seen dynamically as 
we discuss below. 

The sensitivity of the measured SERS spectrum to the probed 
molecular conformation and orientation within the field allows for 
complex and varied molecular flexing and movement to be observed. 
We are able to reliably record SERS trajectories from each individual 
gold NPoM spanning the lipid membrane, for times of order many 
hours, with SERS intensities remaining strong (Fig.S6). In Figure 3a, a 
section of such a time series is seen which shows that far from being 
static over the course of measurement, specific Raman vibrations of 
the lipids undergo fluctuations in intensity and spectral position. 
Analysis of these data sets reveal some clear patterns, with two features 
implying that only few molecules are being probed: (a) opposite sim- 
ultaneous shifts in vibrational frequency of different lines can be seen, 
and (b) a number of lines can appear and disappear simultaneously. 

The Raman intensity for each particular line varies, but in a man- 
ner uncorrelated to other Raman lines. Such observations are con- 
sistent with single-molecule measurements^^'^^ and, although a 
'blinking' mechanism is not yet clear, it may arise from perturbations 
of the molecular excited state lifetime^ ^ or the Raman vibronic coup- 
lings^. More interesting are the simultaneous correlated shifts in 
Stokes frequencies of different Raman modes, seen for the first time. 
The variation between different vibrational lines in the frequency 
shifts and directions strongly implies that vibrations seen here are 
localized to particular parts of each molecule. 

A summary of our measured rms mode shifts on the narrowest 
lines averaged over 1000 spectra (Table 1) reveals shifts of magnitude 



2-6 cm"\ Compression from metal tips in TERS has been suggested 
to shift modes, with estimates that in adenine molecules 10% bond 
contraction yields a 5 cm"^ shift^^'^^'^^. This flexing of the lipid mole- 
cules within the probed gap thus causes both amplitude and fre- 
quency fluctuations. Conventional measurements average over all 
such dynamical and local information, but single-molecule SERS 
using gap-plasmons can reveal these synchronous shifts, further dis- 
cussed below. The ever-present 1380 cm"^ CH2 deformation seen is 
common to all membrane molecules and since these modes are very 
localized in the chains it is least sensitive to larger structural defor- 
mations - we thus use it as normalization for the total number of 
molecules probed. 

We discuss now typical key events (Fig. 3) from many among the 
full data sets (for others see Fig.S7). One dramatic event (typically 
every few minutes) is a sudden change in the entire Raman finger- 
print (Fig. 3a). Assigning the molecular bonds for each peak allows 
this event to be understood as a change in conformation of the lipid 
molecule. Beforehand, clear Raman signatures of the lipid tail, v(C-C 
trans) at 1139 cm"^ and v(C=C) around 1550 cm"\ are seen as well 
as the characteristic mode of the ester linking group v(-[C=0]-0-C) 
at 1183 cm-\ 

Initially the small frequency shifts in these lines are correlated (Fig. 
3a, right), implying they originate from a single lipid molecule. 
During the rapid (<ls) dramatic perturbation of the SERS spectra, 
bands disappear while new bands emerge. These include the char- 
acteristic mode of the phosphate group v(P02) around 800 cm" ^ and 
a clearly shifted trans v(C-C) at 1 120 cm~\ The stretch of the double 
bond v(C=C) centred around 1550 cm"^ persists, suggesting a 
change in the orientation of the lipid head group with respect to 
the optical field polarization in the gap. This change in conformation 
can be followed only because of the precisely defined field orientation 
in our gaps and the well-defined Raman dipole moments of the 
lipids. Chemical transformations may also induce such reconfigura- 
tions, but are not considered likely in the environment here. Along 
with the spectral rearrangement, a strong persistent increase in SERS 
amplitude is also seen which may arise from the improved alignment 
of Raman dipoles with the gap field. 
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Figure 3 | Fast dynamics of hybrid lipid-bilayer (POPCiDOTAP) SERS. (a), Time dependent evolution of SERS from a single gold NPoM. Trajectories 
of individual lines show dramatic jumps and correlated shifts, corresponding to lipid dynamics within the gap plasmon volume. Right: Selected short- 
time dynamics highlighting the (anti-) /correlations between Raman lines, (b), Raman trajectory, including flexing events labelled I, II extracted to 
right, (c). Schematic lipid flexing within gap plasmon hot spot, (d), i, Frequency shifts at each time step of Raman line pairs (vi, V2) showing positive 
(orange) and negative (purple) correlations, ii, Strength of correlation between pairs of Raman lines (vi, V2) with size of points giving correlation strength 
R^. iii, Magnitude of correlation (as in ii) when separation between bonds involved increases (dashed linear fit has R^>0.5). All spectra acquired at 
}^R=633 nm with acquisition time of 1 s. 



More intricate changes in lipid conformation can be observed on 
the timescale of seconds (Fig. 3b). We refer to correlated changes in 
vibrational modes as 'molecular flexing', and explore the changes 
within events I and II. Within I, three bands are seen that exhibit 
synchronous frequency shifts with a simflar magnitude of 10 cm"^ 
(from the previous 1183, 1550 bands and also one at 1070 cm"^), 
suggesting they all report different aspects of a global molecular 
deformation. The line at 1070 cm"^ reports a gauche to trans con- 



formational change of the alkane tail (C-C)^^. Such discrete switching 
is possible only if a single molecule is being interrogated, unless 
several lipids move completely in tandem which is very unlikely 
given the fluidity and activity of the membrane. Changes in leg 
conformation increase the strain of this lipid inside the fluid leaflet 
and this strain is reflected in the downshift of the 1183 cm"^ ester 
linking group. However while this molecular flex takes place, some 
modes of the lipid remain unaffected, such as the 1380 cm"^ CH2 
deformation of the alkane legs. On the other hand, correlated shifts 
for choline, ester, and alkene bonds suggest strain between the head 
group (which is strongly attracted to the nanoparticle) and the tail 
which is 'jittering' in the membrane beneath. 

The SERS trajectory then shows a rapid evolution to fingerprint II 
(Fig. 3b) which is coupled to the events within I, and interpreted as 
arising from the same molecule. Enhancement of the 1325 cm"^ 
methyl deformation (lipid head group) mode likely refiects its 
realignment to the optical field direction, but its simultaneous spec- 
tral shifts are anti- correlated with the carbonyl group stretch 
[v(C=0), around 1700 cm"^]. We find in II clear anti- correlations 
(with spectral line shifts mirrored in opposite directions) in the head 
deformation band at 1320 cm"^ and carbonyl stretch at 1700 cm"^ 



Table 1 | Average flexing mode-shifts for selected assigned Raman 
modes. The central mode frequency and the rms magnitude of 
temporal frequency shifts is measured from a sample of 1 000 
measurements (v = stretch, 6 = bend) 

Mode assignment Centre Frequency (cm"^) rms mode-shift (cm"^) 

5(CH2) 1339 6 

v(-C=0) 1700 6 

v(-[C=0]-0-C-) 1183 2.3 

v(P02) 1239 3.0 
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(Fig. 3d)^^. These spectral dynamics, which can only be explained as 
arising from individual lipids, highlight our capability to now mea- 
sure molecular flexing in real time, and indicate the large number of 
ways strain can be partitioned between parts of a molecule. Our 
ability to collect long time sequences of molecular flexing requires 
advances in theoretical tools to directly extract conformational 
dynamics. However the features highlighted here already clearly sig- 
nal the wealth of information on the membrane that can be adduced. 
We analyse correlations in the spectral shifts between pairs of vibra- 
tions (vi, V2), restricted within shorter (30s) time spans on a single 
molecule, showing in Fig. 3(d,i) both positive (orange) and negative 
(purple) dependencies. We extract all of these correlations given by 
slopes gab = Av^/v^, and capture in Fig. 3(d,ii) this 'correlation fin- 
gerprint' as Cab = sin(2 atan gab) which records how close the two 
frequency shifts are perfectly locked ( + 1) or inverted (-1) or uncor- 
related (0). The size of each Cab point records the correlation coef- 
flcient of the linear regression (how good their correlation is). In 
general we flnd positive correlations are stronger than negative ones, 
reflecting the local response of neighbouring bonds. This can be 
quantitatively addressed by plotting the correlation strength for the 
assigned bonds vs their separation, in terms of number of atoms, in 
Fig. 3(d,iii). Typically we find decreasing correlation strength at more 
remote separations. These correlations extracted from the flexing 
dynamics thus encode detailed information on how molecules move. 

Besides key questions about the rate of molecular flexing (prev- 
iously hidden) for which we are developing faster spectroscopies, and 
how it varies in different lipid layers, fundamental puzzles remain 
about the success of the technique. Brownian motion would be 
expected to move molecules through the probe volume within 
10 |is, hence pinning of the charged DOTAP rafts^^ and interactions 
with the Au nanoparticle must retard lipid motion. Such pinning is 
consistent with the observation of reduced diffusion of membrane 
lipids bearing adsorbed charged polymers on the membrane sur- 
face^^ In addition, similar to recent TERS^^'^^'^\ the lateral resolution 
observed exceeds predictions based on simple models of the optical 
fleld in nanogaps between smooth surfaces and sharply bounded 
metals (Fig.SS). The influence of atomically- sharp corners on the 
facetted nanoparticles may play a role in further fleld localization, 
leading to single nanoparticle SERS over a few lipid molecules. Such 
flelds can also exert optical forces at the nanoscale which influence 
molecular diffusion. 

Vibrational spectroscopy of single molecules within lipid bio- 
membranes presents a powerful means to explore and characterize 
molecular interactions at the nanoscale. Real-time spectroscopies 
will enable greater understanding of the precise interaction of pro- 
teins and lipids in the dynamic and complex cell membrane surfaces. 
Our flndings here open the possibility for monitoring of other cell 
processes relevant to fundamental cellular operation at the single 
molecule level in biological environments. 

Methods 

Sample Preparation. NPoM - HBL formation: A flow cell chamber was fabricated 
from PDMS with a thin 10 [im glass coverslip to provide a viewing window. Gold of 
thickness 50 nm was evaporated onto 10 nm Ti/350 nm SiO2/500 [im silicon by e- 
beam evaporation. The rms roughness was determined to be less than 0.5 nm. The 
ODT monolayer was formed by submerging the substrates in a 5 mM ODT solution 
in dichloromethane (DCM) for 24 h at room temperature. The samples were 
thoroughly rinsed with DCM and dried under nitrogen flow. The contact angles of 
deionized water on this surface (advancing: 110.0° and receding: 104.2°) were 
obtained by a FTA4000 Microdrop goniometer (First Ten Angstroms, Inc). The 
thickness of the ODT SAM (2.1 ± 0.1 nm) on Au was characterized by a UVISEL 
Spectroscopic Ellipsometer (Horiba, Ltd). The rinsed ODT-Au SAM was placed 
within a home-built flow cell and lipid solution flowed into it and left for 1 hour. The 
excess was rinsed with a tris(hydroxymethyl)aminomethane (TRIS) buffer solution 
and the AuNPs (80 nm diameter, BBI Solutions, UK) similarly flowed through, 
rinsed and removed after 1 hour. Within several minutes, accumulation of AuNPs on 
the surface was observed, and the times were adjusted to ensure a suitable surface 
density for the experiments. 



Experimental. Single particle scattering spectroscopy: White light scattering 
spectroscopy was performed in situ on the membrane with fibre-coupled 
illumination from a 40 W halogen lamp (—60° incidence angle). The scattered light 
was collected with a long working distance x50 objective (of the Raman microscope) 
and fibre-coupled into a TEC-cooled OceanOptics QE65000 spectrometer. The 
integration time was 1 s. 

Raman spectroscopy: A Renishaw In Via Raman microscope was used for all 
measurements in point scan mode on the flow cell sample. A x50 long working 
objective was used to collect all spectra. A grating of 1200 lines/mm was used with an 
appropriate edge filter. Excitation of Raman modes was with 633 and 785 nm HeNe 
and diode lasers, typically ~ 1 mW. 
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